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Abstract 

A  systematic  evaluation  of  vertical  velocities  associated  with  stratified  atmospheric  flows 
over  non-homogeneous  terrauis  apparently  has  not  been  reported  in  the  literature.  In  this  study 
we  approach  the  problem  numerically  and  analytically.  Through  a  non-linear  me  del  we  evaluate 
tlie  range  of  the  parameters  for  validity  of  the  linear  theory.  Through  simple  analytical  theory 
wr  estimate  the  role  played  by  the  relevant  parameters. 

Results  Indicate  that,  when  the  transition  in  surface  roughness  is  gradual  between  a  smooth 
and  a  rough  surface,  the  perturbation  of  the  vertical  velocity  heis  the  same  horizontal  scale  as 
the  perturbing  source.  The  nature  of  the  perturbation  depends  on  the  product  between  the 
horizontal  scale  of  the  rough  patch  and  the  Scorer  parameter  of  the  ambient  atmosphere;  i.c. 
very  small  values  of  this  product  (weakly  stratified  atmosphere)  gives  a  wave  trapped  around 
the  top  of  the  stress  layer,  while  values  of  the  order  of  unity  give  a  non-hydrostatic  gravity 
wave  which  propagates  away  from  the  top  of  the  stress  layer.  Values  larger  than  unity  yeld  a 
propagating  hydrostatic  wave. 
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*  W  hen  the  transition  in  surface  roughness  is  abrupt,  the  wave  is  nonhydrostatic  with  a  hor¬ 
izontal  and  vertical  wave  number  equal  to  the  ambient  Scorer  parameter.  When  the  product 
between  the  horizontal  width  of  the  rough  patch  and  the  Scorer  parameter  is  large,  there  is 
no  significant  interference  between  the  wave  upstream  and  the  wave  downstream  of  the  rough 
patch;  the  two  waves  are  of  the  same  nature  and  aUiplitude,  but  opposite  sign.  W'hen  this  prod¬ 
uct  is  small,  however,  the  wave  upstream  interferes  destructively  with  the  wave  downstrear  , 
i.e.  the  amplitude  of  the  resulting  wave  is  accordingly  reduced. 

W’hen  the  surface  roughness  is  periodic,  resonant  amplification  occurs  when  the  horizontal 
wave  number  of  the  perturbation  approaches  the  ambient  Scorer  parameter. 

Since  the  depth  of  updraft  well  exceeds  the  height  of  the  internal  boundary  layer,  this  process 
can  be  important  in  triggering  cumulus  clouds  and  may  have  an  impact  on  the  dispersion  of 
pollutants. 


1  Introduction 

an  air  mass  approaches  a  region  where  there  is  a  substantial  increase  in  surface  roughness, 
v.f'  c>.\i)cct  a  decrease  of  the  air  speed  in  the  lower  layer  and  the  development  of  ascent  associated 
V,  ith  the  resultant  horizontal  convergence.  Such  situations  are  typical,  for  example,  in  coastal 
urban  areas  when  onshore  flow  occurs.  However,  it  may  also  be  of  significance  in  inland  urban 
areas,  because  of  the  contrast  with  surrounding  agricultural  rural  areas  or  when  there  is  a  contrast 
between  prairie  and  wooded  areas.  The  developed  vertical  velocities  may  trigger,  under  supportive 
.synoptic  conditions,  convective  clouds.  The  features  of  the  induced  vertical  velocity  may  also  have 
an  importance  in  dispersing  pollutants. 

In  this  paper  we  approach  the  problem  of  the  vertical  velocity  which  arises  because  of  horizont  al 
inbomogenities  in  the  surface  shearing  stress  in  the  atmospheric  planetary  boundary  layer.  This 
study  is  an  extension  in  more  general  terms  of  a  previous  paper  (Dalu  et  al.  1988),  where  we  report 
on  the  waves  generated  by  change  in  surface  roughness. 

.N  umerous  studies  have  been  carried  out  evaluating  the  impact  of  a  sudden  change  in  the  surface 
roughness  on  horizontal  flow  features  in  a  neutrally  stratified  boundary  layer,  focusing  mainly  on 
the  modification  of  the  flow  w’ithin  the  surface  layer.  In  addition,  considerable  attention  has  been 
given  to  evaluate  th*  growth  of  thermal  and  aerodynamics  internal  boundary  layers.  Hunt  and 
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Simpson  (1982)  provide  an  excellent  review  of  the  studies  reported  by  that  time.  Additional  studies 
reporting  on  the  impact  of  a  sudden  change  in  the  surface  roughness  on  the  horizontal  flow  are 
given  by  Pedergrass  and  Arya  (1984). 

''vertiU,  however,  no  systematic  attention  has  been  given  in  the  evaluation  of  the  features  of 
the  vertical  velocities  generated  tn  a  stable  stratified  flow  when  the  air  crosses  a  change  in  the 
sur.'^ace  roughness  (CSR).  Wagner  (1966)  reported  on  the  structure  of  the  vertical  velocity  induced 
by  a  CSR  as  computed  through  a  numerical  model.  However,  he  neglected  to  discuss  the  vertical 
e.xtensions  of  the  updraft  and  its  wave  characteristics,  because  of  his  upper  boundary  condition:  in 
fact  he  required  the  vertical  velocity  to  be  zero  and  the  flow  to  be  in  geostrophic  balance  at  1000 
meters  above  the  ground. 

Claussen  (1986)  computed,  using  a  model  simulation,  the  vertical  velocity  due  to  a  CSR.  How¬ 
ever,  the  computed  vertical  velocity  was,  in  general,  very  sensitive  to  the  horizontal  grid  resolution, 
whicli  should  be  reduced,  in  many  occasions,  to  several  hundred  meters  in  order  to  appropriately 
resolve  the  related  vertical  velocity. 

Using  a  very  coarse  horizontal  resolution,  Vukovich  and  Dunn  (1978)  in  their  numerical  model 
simulation  of  the  St.  Louis  urban  area,  suggested  that  the  surface  roughness  has  only  a  small  effect 
on  the  circulation  for  the  wind  speeds  used  in  their  study.  Alestalo  et  al.  (1985),  using  a  hydro  tatk 
two-dimensional  model  with  a  grid  interval  of  4  km  simulated  the  airflow  in  the  Baltic  shore  region  of 
Finland,  and  found  a  maximum  for  the  vertical  velocity  of  order  of  1  cm  ,  due  to  the  CSR.  They 
attributed  the  reported  increase  of  precipitation  in  that  area,  in  the  absence  of  thermal  forcing,  to 
ti.*'  vertical  velocity  induced  by  the  CSR.  Pielke  (1974)  evaluated  the  magnitude  of  vertical  velocity 
cau.scd  by  a  CSR  over  Florida  using  all  km  horizontal  resolution  model.  Although  the  magnitudes 
were  small  0.1  cm  s~*),  it  was  concluded  that  shallow  warm  rain  cloud  over  the  southeast 
coast  of  Florida  could  be  due  to  this  mechanism.  Finally,  Roeloffzen  et  al.  (1986)  presented  a 
steady  state  model  calculation  of  secondary  flow  patterns  forced  by  a  CSR.  Adopting  a  neutral 
boundary  layer  and  using  a  refined  grid  resolution,  they  suggest  that  frictional  effects  involved  with 
a  CSR  at  a  coast  line,  can  lead  to  a  secondary  circulation  on  the  mesoscale.  They  suggest  that  this 
forcing  is  a  factor  in  the  observed  coastal  frontogenesis  active  in  the  early  fall  along  the  coast  of 
the  .N'etheriands. 
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The  present  study  presents  a  systematic  analytical  and  numerical  evaluation  of  CSR  impact 
in  a  stably  stratified  air  mass,  for  an  one-layered  atmosphere  and  in  a  two-layered  atmosphere, 
on  the  generation  of  vertical  velocities  and  their  characteristics  as  dependent  on  key  parameters 
including  the  wind  speed,  the  thermal  stratification  of  the  lower  atmosphere,  and  the  magnitude 
of  the  shearing  stress. 

The  results  of  this  study  are  confirmed  by  the  general  mesoscale  features  observed  during  the 
-Amazon  Boundary  Layer  Experiment  ABLE  2B  (Garstang  et  al,  1990).  However,  the  present 
theory  should  be  extented  to  include  the  diabatic  sources,  in  order  to  explain  the  ABLE  2B  results 
in  a  more  quantitive  way. 

2  Waves  Generated  by  a  CSR  as  Simulated  by  a  Numerical  Model 

Here  we  present  some  results  concerning  the  vertical  velocity  induced  by  a  sudden  change  in  terrain 
surface  roughness,  as  simulated  by  a  numerical  model.  For  the  simulation  w'e  use  the  non-hydrostatic 
version  on  an  /-plane  of  the  CSU  Regional  Atmospheric  Meteorological  Model  RAMS  (Tremback 
cm  et  al,  19....???).  The  ambient  flow  is  perpendicular  to  the  CSR,  and  the  atmosphere  is  stable 
stratified.  We  assume  a  sudden  change  of  roughness  Zo  =  I  rn,  and  a  flow  intensity  of  1/  =  5  m 
and  U  =  10  m  s"*.  We  show  the  results  after  a  few  hours  of  simulation,  i.e.  when  the  system 
evolves  to  an  almost  steady  state.  Results  for  an  ambient  flow  of  V  —  b  m  s"*  are  shown  in 
Fig.i  and  for  an  ambient  flow  of  U  =  10  m  in  Fig.2.  The  results  in  Fig.l  and  in  Fig.2  show 
that  substantial  vertical  velocities  can  be  induced  by  a  sudden  change  of  surface  roughness  and 
that  the  perturbation  is  of  a  wave  nature  (since  the  atmosphere  is  stratified)  with  a  wave  length 
and  intensity  which  doubles  for  a  doubling  of  the  wind  intensity.  This  response  suggests  a  linear 
behavior. 

The  non-linear  behavior  is  apparently  confined  to  within  the  stress  layer,  which  in  the  numerical 
model  is  parameterized  according  a  turbulent  kinectic  energy  TKE  scheme  (Deadorf,  1980).  As 
a  cosequense,  the  depth  of  the  stress  layer  and  the  shear  velocity  do  not  respond  linearly  to  the 
variations  of  the  flow  intensity. 

However,  from  the  simulations  completed  with  the  nonlinear  model,  we  found  that  for  AUfU  < 
.2  with  the  largest  perturbation  near  the  surface,  i.e.  the  relative  perturbation  of  the  horizontal 
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velocity  is  less  than  20%  (similar  results  have  been  found  by  Claussen,  1987).  Thus  we  deduce 
that  the  response  in  the  atmosphere  above  the  stiess  layer  is  linear  and  that  the  non-linearities  are 
confined  to  the  turbulence  within  the  stress  layer. 

We  proceed  in  investigating  the  nature  of  the  perturbation  through  linear  analytical  theory,  in 
order  to  find  simple  relations  between  the  behavior  shown  in  Fig.l  and  in  Fig. 2  and  the  ambient 
parameters.  In  the  following  theoretical  study,  we  assume  the  depth  and  the  intensity  of  the  stress, 
whicli  ran  be  obtained  through  observations  or  turbulence  theory. 


3  The  Governing  Equation  for  the  Linear  Problem 

If  we  assume  that  the  process  is  stationary,  two-dimensional  and  Boussinesq,  then  the  primitive 
equations  in  linear  form  can  be  reduced  to  a  Scorer  type  equation  for  the  vertical  velocity  in 


non- homogeneous  form; 


The  linearization  is  justified  by  fact  that  the  non-linear  numerical  simulations  have  shown  that  the 
perturbations  are  small  in  comparison  to  the  intensity  of  the  ambient  wind,  see  section  (2).  For 
a  derivation  of  equation  (3.1),  see  Eliassen  (1977).  The  hat  denotes  the  Fourier  transform  of  the 
variable,  k  is  the  horizontal  wave  number,  kg  =  f /U  is  the  inertial  wave  number,  (/  i-  the  Corioli:^ 
parameter,  U  is  the  ambient  flow  perpendicular  to  the  change  in  surface  roughness,  and  bz  is  it' 
shear),  r  is  the  resulting  shear  stress,  u'  is  the  shear  velocity,  and  I  is  the  Scorer  parameter  (Scorer. 


1953): 


,2  .  u  ,.r2 

'  =z?T-ir 


.vhere  N  is  the  Brunt-Vaisala  frequency  and  h  is  the  buoyancy  of  the  environmer..  Equation  (3.1) 


can  rewritten  as; 


-b  u\k)<i  =  G^k)^  =  G\k)  ^ 


with  i/2(A:)  =  ^  and  ^  p-^Tp 


5 


In  the  wave  number  region  where  i/^{k)  <  0,  the  waves  are  trapped  around  the  perturbing  source 
within  an  e— folding  vertical  distance  equal  to  /<”*•  The  vertical  wave  number  Ho  for  the  trapped 


waves  is: 


Ho{k)  =  |i  i’(k)j  =  wncn  u  |ft.|  v  «■<>  vi  wucn  •  v,  |n,| 

In  the  wave  number  legion  where  v^ik)  >  0,  the  w’aves  propagate  away  from  the  perturbing  source 
with  a  vertical  wave  number  equal  to  fii.  The  vertical  wave  number  n\  for  the  propagating  waves 


1/2  _  k-i 


when  0  <  |i(:|  <  ko  or  when  /  <  (A:j  <  oo 


//2  _  1-2 

fiiik)  =  t/(k)  =  kJ when  ko  <  |*:|  <  /. 


3.1  The  Stress  and  the  Shear  Velocity 

For  simplicity,  we  assume  that  the  stress  has  the  same  direction  and  opposes  the  ambient  flow. 
Furthermore,  since  the  stress  is  assumed  to  decay  linearly  with  altitude  within  the  stress  layer  (this 
simple  choice  is  suggested  by  the  results  of  the  numerical  simulations;  i.e.  see  section  (2)): 

t(x,2)  =  To  ne{h  —  z) — ^ — F{x)  (3T)) 

where  Tq  is  the  surface  shear  stress,  F{x)  is  its  horizontal  distribution,  and  He  is  the  Heaviside 
function.  We  study  the  atmospheric  response 
to  a  jicriodic  horizontal  distribution  of  the  stress: 

case  (i)  cos(q3:)  =»  fz^{k,z)  =  - 


to  a  bell  shape  horizontal  distribution  of  the  stress: 


,  ToH^-h)  X  Tof{z-h)7ra 

case  (ii)  r,j(i,z)=  - — - ^2^^772  ^ ~  - ^exp(-al),  (3.8) 

and  to  a  top  hat  horizontal  distribution  of  the  stress: 

case  (ill)  r„(i,2)  =  [He{x  +  a)  -  tfe(i  -  a)]  =>  f^^{k,z)  =  - ^—^—^.(3.9) 

Here  f'  is  the  Dirac  function.  The  tilde  denotes  the  cosine  Fourier  transform: 


F(k)  =  f 
Jo 


dxF(x)cos{ki) 
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The  relation  between  the  stress  r,  the  wind  intensity  U ,  the  surface  drag  Co,  and  the  shear 
velocity  u'  is  given  by  (Panofsky  and  Dutton,  1984): 


... 


(3.10) 


Here  k  is  the  von  Karmann  constant  and  Zq  is  the  surface  roughness.  The  order  of  magnitude  of 
the  depth  of  the  stress  layer  is  (Blackadar  and  Tenneks,  1968): 


h 


so  that 


Ih  0.3 


lu' 

T 


0.3  Cp 


/ 


(3.11) 


These  simple  formulas  give  a  shear  velocity  and  a  stress  layer  which  grow  almost  linearly  with  the 
wind  intensity,  with  a  depth  almost  a  constant  fraction  of  the  inverse  of  the  Scorer  parameter. 

However,  results  from  the  non-linear  model  (section  2)  show  that  the  the  shear  velocity  and 
the  stress  layer  grow  almost  linearly  with  the  wind  intensity  when  the  wind  is  weak;  their  growth 
is  slower  than  linear  at  a  higher  wind  speed.  The  quantification  of  this  result  is  rather  sensitive  to 
the  turbulence  scheme  used. 


4  Atmospheric  Response  to  a  Roughness  Changes 


Tlie  atmospheric  response  to  a  general  surface  roughness  change  is  given  by: 


W{l.  Z)  —  -f-  /i  -t-  /oj 


j  dkGo{k)w^Ak,z)Fik)  +  I  dkG,(k)w^,{k,z}r(k) 

0  *o 


-I-  J dkGo{k)w^Jk,z)F(k) I  with  =  'yTh 

Where  u-  is  thf’  amplitude  of  the  perturbation  of  the  vertical  velocity.  Here 

,  _  IG^jk) _ Ik  _  lG‘(k)  _  Ik 

^o{k)  y(krzrkTfii2GrW)  fii(k)  -  kl)(F  -  k^) 

The  u'^^Jk,  z)  waves  are  trapped  around  the  top  of  the  stress  layer: 

U'oik.z)  =  {H€(z  -  h)exp{-fio{:  -  h))  -I-  (^c(r)  -  He{z  -  /i)] e.xp( -//^(/i  -  z)) 


(4.2) 


-^e(z)exp(-/io(z  +  /i))}  cos(Jtz)  (-1.3) 

In  eq.  (4.3)  the  first  and  the  second  term  give  waves  which  are  trapped  over  and  below  the  top  of 
the  stress  layer,  respectively.  Since  these  waves  have  their  maximun  amplitude  at  the  top  of  ihe 
stress  layer,  they  may  have  a  role  in  triggering  cumulus  convection,  even  if  their  amplitude  decays 
exponentially  with  distance;  the  third  term  is  the  wave  reflected  from  the  ground,  (see  Appendix). 
The  w^j{k,z)  waves  propagate  away  from  the  top  of  the  stress  layer: 

=  {He[z  -  /i)sin(pi(2  -  h)  +  kx)  +  [He{z)  -  He{z  -  h)]s\n{^\(h  -  2)  +  kx) 

—He{z)sm{fii{z-\-h)  +  kx)}  (4.4) 

In  eq.  (4.4)  the  first  term  is  the  wave  which  propagates  upward,  the  second  term  is  the  wave  which 
propagates  downward,  and  the  third  term  is  the  wave  reflected  from  the  ground,  (see  Appendix). 

We  assume  the  following  val'  *^''r  environment  parameters  (when  the  atmosphere  is  stratified): 
A  =  10“^  see"*,  U  =  10  m  u"  =  20  cm  and  /  =  10~^  s~’,  then  we  have  /  =  10~^  m~^ , 
k\  =  10~®  m~‘,  h  =  300  m\  and  A'  =  /  =  0.0  when  the  environment  is  non-stratified. 

4.1  Resonant  amplification  due  to  periodicity,  (case  (»)) 

W  hen  the  surface  roughness  is  periodic  (eq.  (3.7))  the  vertical  velocity  is  given  by: 

U’(l,2)  =  tI-{Go(Q)U'^„(Q,2)  +Gl(Q)uv,(a,2)}  (4.-5) 

Where  Go(o)  =  — >  ■■■  .  and  Gi(q)  =  (4.6) 

x/(A-2-a2)(/^-Q2)  ^  ^(a2-^2)(/2-a2) 

with  u'o(q,2)  =  -  {fle(z  -  h)exp(-fXo(z  -  h))  +  [He{z)  -  He(z  -  h)]exp(-fio{h  -  2)) 

-He(z)^xp{-po(z  +  h))}  cos(aa-)  (4.7) 

a!Kl  with  u-^j(a,2)  =  {He{z  -  /i)sin(pi(2  -  h)  +  or)  +  \He{z)  -  He{z  —  h)]sm{pi{h  -  2)  +  qj) 


■He{z)s\n(p^{z  +  h)  +  or)} 
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fio  =  a 


-  q2  ^  //2  -a2 

kl-a^ 

The  amplitude  Go(a)  of  the  trapped  wave  and  the  amplitude  Gj(a)  of  the  propagating  wave  are 
amplified  by  resonance  when  a  =  I,  i.e.  Go{q  —*  1)  -*  oo  and  Gi{a  -+/)—►  oo.  When  the  wave 
number  a  of  the  rough  patch  distribution  approaches  the  ambient  Scorer  parameter  I,  the  maximum 
enhancement  of  the  vertical  velocity  occurs. 

There  is  a  resonance  also  at  q  =  ko,  but  this  resonance  is  canceled  by  infinitely  rapid  oscillations 
Pi(a  ^  ►  00,  or  by  infinitely  strong  trapping,  /io(o  -♦/)—►  oc. 


4.2  Vertical  velocity  excited  by  a  bell  shaped  stress,  (case  (u)) 

in  this  section  we  study  the  vertical  velocity  induced  by  a  rough  patch,  in  relation  to  its  horizontal 
extension.  The  bell  shape  distribution  is  ideal  for  this  kind  of  analysis,  as  shown  by  Queney  (19-17) 
and  by  Smith  (1979)  for  the  vertical  velocity  induced  by  a  bell  shape  mountain.  For  a  bell-shaped 
distributed  stress  case  (ij),  from  eq.  (4.1)  the  vertical  velocity  is  given  by: 

u(i,  z)  =  /<,,  -f  /i  +  /oj  =  (^kGo{k)u-^^{k,  z)  a  exp{-ak) 


J  dkGi{k)tL'^^{k,  z)  a  exp{-ak) j  ° 


Fig. .3a  shows  the  atmospheric  response  to  a  bell  shape  horizontal  distribution  of  the  stress  with 
(1  =  8  km,  {la  =  8).  Fig.4a  shows  the  atmospheric  response  to  a  bell  shape  horizontal  distribution 
of  the  stress  with  a  =  1  km,  {la  =  1).  The  genera!  features,  for  standard  values  of  wind  velocity 
and  atmospheric  stratification,  are  that  the  vertical  wave  number  equals  the  Scorer  parameter  /. 
the  horizontal  extension  of  the  perturbation  is  of  the  order  of  the  width  a,  and  the  propagating 
wave  dominates  the  trapped  wave. 

Let  us  examine  in  more  detail  the  atmospheric  response  for  different  widths  a. 

Some  useful  approximations 

(a)  -  When  kgO  K  1  or  larger,  /q,  >>  7]  and  the  wave  is  trapped 

Due  to  the  exponential  decay  of  the  Fourier  transform  of  the  bell  function  for  increasing  values 
of  the  wave  number,  when  the  rough  patch  is  very  very  large  (i.e.  when  the  horizontal  scale  is 
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comparable  to  the  size  of  the  inertial  wave),  the  contributions  of  the  second  and  third  integrals  are 
negligible  in  comparison  to  the  contribution  of  the  first  integral; 

w{x,z)  ss  /o,  =  ti)  /  dkGo(k)Wft^{k,  z)  a  exp(-ak) 

Jo 


{ 


He(z  -  h) 


(a  +  z  -  h)^  — 

[(a  +  z  - /i)2  +  z2)2 


+  [He{z)  -  ne{z  -  h)] 


{a  +  h  —  —  x"^ 

[(a  +  /i  -  2)2  +  2r2]2 


JL  I  =  k  and  Go(k)  =  -^  (4.10) 

[(a  +  2  +  hy  +  i2]  J  ko 

W'hich  is  a  trapped  wave.  The  vertical  velocity  monotonically  decreases  with  altitude  above  the 
stress  layer,  and  its  structure  does  not  depend  explicitly  on  stabilty  parameters  (although  h  may 
depend  on  them),  as  in  the  Ekman  solution. 

(b)  -  When  la  >>  1  (but  kga  <  IJ,  Ii  >>  +  /©j,  and  the  wave  is  hydrostatic 

The  trapped  wave  contribution  is  negligible  in  comparison  to  the  contribution  of  the  propagating 
hyd;'“tatic  v.avt; 

fOO 

V'(x,z)  s;  /j  =  tZ’  /  dkGi(k)w^^(k,z)  a  exp(-ak) 

Jo 

=  ^  {He{z  -  h)  [asin(A(2  -  h))  +  icos(A(2  -  h))]  +  [He{z)  -  Heiz  -  h)] 


[a  sin(A(/t  -  2))  +  I  cos(A(/i  -  2))]  —  He{z)  (asin(A(2  +  h))  +  z  cos(A(2  +  h))]}  (-l-l  1) 


with  =  /;  G\(k)  =  1  and  A  =  -  kl 

When  the  atmosphere  is  stratified  and  the  rough  patch  is  large  (with  no  abrupt  roughness  tran¬ 
sitions),  the  vertical  velocity  has  a  hydrostatic  wave  structure  with  a  vertical  wave  number  equal 
to  the  Scorer  parameter  (corrected  because  of  inertia)  and  a  horizontal  scale  comparable  with  tlic 
width  of  the  rough  patch.  Fig.3b  shows  the  hydrostatic  waves  for  la  =  8,  (eq.  (4.4)),  and  how  well 
this  form  approximates  the  full  solution,  (eq.  (4.9)),  shown  in  Fig.3a. 

(c)  -  IlTjcn  la  =  0(1),  lo,  <<  I\  -I-  /03,  and  the  wave  is  non-hydrostatic 
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When  the  atmosphere  is  stratified,  but  the  rough  patch  is  small,  the  inertial  effect  is  negligible: 
u'(i,r)  w  /o,  +  /o2  =  tu  dkGi{k)w^^{k,z)  a  exp(— fca) 


dkGoik)w^^(k,z)  a  exp(-/:a)| 


(4.12) 


The  propagating  wave  is  non-hydrostatic,  with  a  wave  number  equal  to  the  Scorer  parameter: 

/i  =  u'  la  exp(-/a)-  {IIe(z  -  h)[sin{l{z  -  h))  Jo{l{x  -  z  +  h))  +  cos{l{z  -  h))  Ho{l{x  -  r  +  h))] 


+  =  )  -  He{z  -  fi)]  [sin(/(/i  -  z))  Jo{l{x  -|-  2  -  fi))  +  cos(/(/i  -  z))  Ho{l{x  -f  2  -  /i))] 


-//€(2)  [sin(/(2  h))  J,{l(x  -  z  -  fi))  +  cos(/(z  +  h))  Ho{l{x  -z-  fi))]} 


with  /ii  =  /  -  A:  and  = 


’ 


Tlie  near  field  contribution  of  the  trapped  wave  comes  from  large  wave  numbers  and  could  be 
written  in  terms  of  the  incomplete  1'  functions,  but  a  good  simple  approximation  is: 

r  -  ,  I  \  {  rr  ,  +  Z  -  h)  r  rr  /  %  r,  ,  .  m  ufc  +  /}  -  z) 

=  «■  exp(-(a)  {«,(,  -  -  «,(,  -  Ml 


[(a  +  fi  -  2)2  -f  j21 


-He{z] 


a{a  +  z  +  h) 
[(0  +  z  +  hy  +  i2 


I  with  Ho  =  k  and  Go(k)  =  1 


The  perturbation  is  confined  in  space,  in  a  region  above  the  rough  patch  at  an  altitude  equal  to 
the  depth  of  the  stress  layer. 

Since  the  far  field  contribution  of  the  trapped  wave  comes  from  wave  numbers  close  to  the 
Scorer  parameter,  the  trapped  wave  contribution  is: 

f  i  =  -y  la  exp(-/a)—  A"o(/x)  {Ht{z  -  h)exp(-l{z  —  h))  -|-  [/fe(2)  -  He{z  -  /i)]exp(-/(h  -  2)) 


■//e(2)exp(-/(2 -t- fi))  with  =  /  and  Go{k)  =  -~==='^ 


11 


Tlie  far  field  contribution  of  the  trapped  wave  decays  exponentially  with  the  distance  from  the  top 
of  the  stress  layer,  therefore  interferes  with  the  propagating  wave  only  at  z  w  /i,  where  the  trapped 
wave  cancels  the  propagating  wave  upstream  and  strengthens  the  propagating  wave  downstream 
(the  Neumann  function  is  even,  while  the  Struve  function  is  odd,  with  the  same  asymptotic  absolute 
value),  producing  a  wake  of  secondary  cells  downstream  at  the  level  of  the  top  of  the  stress  layer. 

Jo,  Ho  and  A'o  are  the  zero  order  Bessel,  Struve  and  Neumann  functions,  respectively,  witli  the 
usual  convention  that  even  function  of  negative  argument  equals  the  function  of  the  absolute  value 
of  the  argument,  and  that  an  odd  function  of  negative  argument  equals  the  negative  of  the  function 
of  the  absolute  value  of  the  argument. 

Summarising,  when  the  atmosphere  is  stratified  and  the  rough  patch  is  small,  the  vertical 
velocity  has  a  wave  structure  with  a  horizontal  and  a  vertical  wave  number  equal  to  the  Score: 
parameter,  with  an  exponeiitially  decreasing  amplitude  for  a  decreasing  width  of  the  rough  patch. 

Fig.4b  shows  the  non-hydrostatic  wave  for  la  =  \  and  how  well  approximates  the  full  solution, 
(eq.  (-1.9)),  shown  in  Fig.4a. 

(d)  -  When  la  <<  1  (but  exp{-koa)  1  >>  /o,  -f  I.,  and  the  flow  is  irrotatioval 

When  the  stratification  is  weak,  /  is  small,  and  the  contribution  of  the  first  and  the  second 
intecrahs  in  eq.  (4.9)  can  be  neglected,  only  the  third  integral  coi.libutes.  Through  the  use  of 
knowledge  that  the  vertical  scale  of  the  perturbation  equals  the  depth  of  stress  layer,  we  simply 

ha  vf': 


,  f  f  ,r  ,  .N  a{a+Z-h)  ,  rrr  r  X  ri  i  ..  M  0(0  +  /j-zl 

lo,  =  —  ^  H e{z  -  h),;-^  ^  jr  +  [^c(z)  - //e(z  -  h)] 


i 


[(a  -h  z  -  h)‘  +  ar‘] 


[{a  +  h  -  zfl  -t  x^] 


/fe(z)  lie  =  k  and  wGo{k)  =  ^ 

'((a-Pz-f /i)2-t-i2]j  ^  '  h 


(4.13) 


The  vertical  velocity  monotonically  decreases  with  altitude  above  the  stress  layer.  Fig. 3c  shows  the 
tr.ifqied  waves  for  /  =  0.0  and  a  =  8  km,  (eq.  (4.13)).  The  trapped  wave  (non-stratified  atmosphere) 
sliiiwn  in  Fig. 3c  considerably  differs  from  the  propagating  wave  (stratified  atmosphere)  shown  in 
Fiy.-ia.b. 
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4.3  Vertical  velocity  excited  by  a  top  hat  stress,  case  (tn) 

In  this  section  we  study  the  vertical  velocity  induced  by  a  rough  patch  of  finite  horizontal  extension, 
with  very  sharp  transitions  between  the  smooth  and  the  rough  region  and  viceversa  ({iii)  in  eq. 
3.9),  this  is  equivalent  to  problem  solved  by  Lira  (1943)  for  the  flow  over  a  step  function  mountain. 
From  eq.  (4.1)  the  vertical  velocity  is  given  by: 

wix,z)  =  ^  I  dkFoik)Goik)  +  J‘  dkFi{k)Giik)  +  dA:F<,(A:)G',(l:)|  (4.14) 


where 


^  sin(A-a)  ^  sin(A:(i  +  a))-sin(A:ljr- a))  , 

^o\k)  =  — - — w^^{k,z)  =  - — - {He{z  -  h)exp{~Ho{z  -  h)) 


+  [He{z)  -  He(z  -  /i)]exp(-/io(h  -  z))  -  H€[z)exp(-fioi^  +  /j))} 


(4.15) 


fiik)  =  w^^{k,z)  =  ^  {He{z  -  /i)[cos(/Uj(2  -  /i)  +  k{x  -  a))  -  cos(fii{z  -  h)  +  k{x  +  c))] 

+  [Heiz)  -  He{z  -  /i)]  [cos(/ii(/i  -  z)  +  k{x  -  o))  -  cos(//i(/i  -  i)  +  k{x  +  a))] 


-He{z)[cos(ni{z  +  h)  +  k{x  -  a))  -  cos(/ii(2  +  h)  +  k(x  +  a))]} 


(4.16) 


From  eq.  (4.15)  and  eq.  (4.16)  we  see  that  the  wave  at  the  upstream  transition  and  at  the 
downstream  transition  are  of  the  same  nature,  but  with  an  opposite  sign.  Furthermore,  when 
'21a  =  0(1),  the  wave  at  the  upstream  transition  (from  smooth  to  rough,  x  =  -a)  interferes 
destructively  with  the  wave  at  the  downstream  transition  (from  rough  to  smooth,  x  =  a);  and 
when  2/a  <<  1  the  perburbation  becomes  negligible,  because  destructive  interference. 

The  resulting  vertical  velocity,  above  the  abrupt  transition  between  the  smooth  and  the  rough 
surface  (Fig. 5a),  is  given  by: 


w(x,z)  =  /o,  +  /i  +  /o- 


2i7'  f  /■*» 

=  —\  dkF;,{k)Go{k)+  /  dkFl(k)G,(k) 

tr  [Jo  Jka 

1"  dkF:,{k)Go{k)^ 
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where 


Ki^')  =  -  h)exp{-^lo{z  -  h))  +  [ne{z)  -  He{z  -  h)]exp{-^lo{h  -  z)) 

-Heiz)exp{~fioiz  +  h))}  (4.18) 

^i(^')  =  -  h)cos{(ii{z  -  h)  +  k{x  +  a))  +  [He{z)  -  He{z  -  /i)] cos(/ii(/i  -  2)  +  k{x  +  a)) 


-He{z)cos(fii{z  +  h)  +  k(x  +  a))}  (4.19) 

From  Fig. 6a,  we  deduce  that,  when  the  transition  between  the  smooth  and  the  rough  surface 
is  abrupt,  the  wave  is  non-hydrostatic  with  a  horizontal  and  a  vertical  wave  number  equal  to  the 
Scorer  parameter.  Since  the  horizontal  scale  of  the  inertial  wave  is  much  larger  than  the  inverse  of 
the  Scorer  parameter,  inertia  is  negligible. 

Some  useful  approximations 

foj  Have  in  a  stratified  atmosphere,  /<,,  <<  Ii  -f 

The  contribution  of  the  first  integral  in  eq.  (4.17)  is  negligible: 

u-ix,z)  s:  -f  7,^  =  ^  y‘ dkF{ik)Gi{k)  +  dkF;,ik)Go{k)^  =  F  +  F  (4.20) 

Til'  propagating  wave  is  non-hydrostatic,  with  structure  similar  to  a  lee  mountain  wave: 

7,  =  I  {77c(2  -  /i)  [sin(/(2  -  h))  HoiKx  +  a  -  z  +  h))  -  cos(/(2  -  h))  Jo(l{x  +  a  -  2  -|-  h))] 

+  [He{z)  -  I{e{z  -  h)]  (sin(/(/i  -  2))  Ho{l(x  +  a  +  z  -  h))  -  cos(/(/i  -  2))  Jo{l{x  +  a  4  z  -  h))] 


-IIe{z)  [sin(/(2  4-  h))  Ho{l{x  a  -  z  -  h))  -  cos{l{z  4-  h))  Jo(l{x  4-  a  -  2  -  h))]} 


Pi  =  I  —  k  and 


Giik)  ^  1 

k  y/P  -  k^  ’ 


The  near  field  contribution  of  the  trapped  wave  is: 


—  (/7c(2  - /i)  tan  ^  -p  [He{z)  -  He(z  -  h)]  tan  * — 

Trf  z  —  h  a  —  z 


/7c(2)tan  ^ 


a-  -I-  a  ) 

TThJ 
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with  f^o  =  k  and  Go(k)  =  1 


The  trapped  wave  contribution  is  confined  above  the  smooth-rough  transition  zone,  witliin  an  h 
distance  from  the  top  of  the  stress  layer. 

In  the  far  field,  the  trapped  wave  contribution  is: 

w 

=  sign{x-\-a)  JoUir  +  a)){He{z  -  h)exp{-l{z  -  h)) 

+  [He{z)  -  He{z  -  /i))]exp(-/(/i  -  z))  -  H e{z)  exp{-l{z  -t-  h))) 

Since  tlie  far  field  contribution  comes  from  wave  numbers  close  to  the  Scorer  parameter,  we  have 


assumed; 


=  /  and 


Goik) 


-  /2 


The  trapped  wave,  in  the  far  field,  decays  exponentially  with  the  distance  from  the  top  of  the 
stress  layer,  therefore  interferes  with  the  propagating  wave  only  when  z  ss  h,  where  the  trapped 
wave  cancels  the  propagating  wave  upstream  and  strengthens  the  propagating  wave  downstream, 
producing  a  wake  of  secondary  cells  downstream  at  the  level  of  the  top  of  the  stress  layer. 

Fig.Sb  shows  how  well  (4.20)  approximates  the  upstream  transition  shown  in  Fig. 5a. 
f’l)  Trapped  wave  in  a  non-stratified  atmosphere,  lo^  >>  hi  +  h 

When  the  atmosphere  is  non-stratified,  most  of  the  contribution  comes  from  the  third  intcrgral 
in  eq.  (4.17)  and  a  good  approximation  of  at  the  abrupt  transition  between  the  smooth  and  the 
rough  surface  is: 

,2  . 

u(x,  z)  s;  /oj  =  (  He{z  -  h)  tan“'  £jFa  ^  _  }Je(z  -  h)]  tan~’ 

TTU  \  Z  —  h  il  —■  z 


-j-  1 

-/fe(z)tan“* - >  with  Po  =  k  and  wGo(k)  = —— 

z  +  h  )  U 


(4.21) 


Fig. 5c.  shows  the  trapped  wave. 
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4.4  Waves  in  a  two  layers  atmosphere 

The  wave  emerging  from  the  stress  layer  is  matched  with  the  wave  in  the  free  atmosphere  using 
the  boundary  value  Green  function  defined  in  the  Appendix. 

(a)  iraue  in  the  free  atmosphere,  generated  in  a  non-stratified  stress  layer 

When  the  stress  layer  is  non-stratified,  I'  =  0,  and  /  ^  0  in  the  free  atmosphere,  the  vertical 
velocity  is: 

f  h\s\n(l{x  -ha))  -  sin(/(^  -  h))]  -  (x  +  a  -  2  fe)[(l  -f  h/)cos(/(j  +  a))  -  cos(/(z  -  h))] 

\  {x  +  a  -  z  -i-  h}'^ 


+h(x  -I-  a) 


(i  +  a)^  +  (2  -  —  2/i(r  -  h) 

[(i  +  a)2-K^-h)2f 


I 


(4.22) 


Fig.O  shows  the  trasmitted  wave. 

fb)  I! 'm  e  ducted  u'ithin  the  surface  layer 

When  tiie  stress  layer  has  a  Scorer  parameter  with  /'>>/%  0,  where  /  is  the  Scorer 
parameter  of  the  free  atmosphere,  the  wave  is  ducted  between  the  ground  and  the  top  of  the  stress 
layer.  The  vertical  velocity  is: 

w{x.  r)  w  —  (  He{z  -  /i)tan~'  ^  [He{z)  -  He{z  -  h)]tan~'  ^4^  _  1 

tI.  z  —  h  h  —  z  z  +  n ) 


-t-^^e(2  -  h)  exp(-/(z  -  h))  Y,  (sin(2m/'/i)  Ho{l'{x  +  a  -  2mh)) 

^  *n=0 


-  cos{2ml'h)  Jo[l'(x  -fa  -  2m/i))  —  sin(2(m  4-  l)l'h)  Ho{l'{x  4-  a  -  2(7n  4-  l)h) 


4-  cos(2(m  4-  l)l'h)  Joil'ix  +  a  —  2(m  4-  l)h)] 
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(4.23) 


+  ^  [H€{z)  -  Ht{z  -  /i)]  ^  [siii(/'((2m  +  l)/i  -  z))  +  a  +  z  -  (2m  +  l)/i)) 

-  cos(/'((2m  +  l)/i  -  z))  Jo(l\x  +  a  +  z  -  (2m  +  l)/i)) 

-  cos(/'((2m  +  l)h  -  z))  Jo{l'{x  +  a  +  z  —  (2m  +  l)/i)) 

+  cos(/'(2  +  (2m  +  l)h))  Jo(l'ix  +  a  -  z  -  (2m  +  l)/i))] 

The  trapped  wave  contributes  only  near  the  smooth-rough  transition,  because  its  amplitude  de¬ 
cays  very  rapidly  with  the  number  t.  reflections  m.  The  wave  is  fully  reflected  from  the  ground, 
penetrates  in  the  non-stratified  free  atmosphere  a  distance,  and  then  bounces  back  towards 
the  ground.  Summarising  when  the  stratification  in  the  stress  layer  is  much  larger  than  the  one  in 
the  above  free  atmosphere,  the  wave  is  ducted  between  the  ground  and  the  top  of  the  stress  layer. 
Fig. 7  shows  the  wave  ducted  within  the  stress  layer. 


5  Subgrid-Scale  Parameterization 

In  section  4.1  we  have  shown  that  resonant  amplification  can  occur  when  the  distribution  of  tlie 
surface  stress  is  periodic.  However  it  is  more  realistic  to  assume  that  the  flow  encounters  a  sudden 
change  of  surface  roughness,  followed  by  smooth  patches  alternated  with  rough  patches: 

_  —  ft)  >  I  .\n  II  ^  /  /  r  ^  \ 


r,i(x,z)  = 


J3(-l)”/fe(j  -  71  a) 


Fig. 8  shows  the  resulting  vertical  velocity  for  differen'  value  of  the  product  la,  obtained  using  eq. 
(4.20).  Results  shown  in  Fig.8  contradicts  the  resonance  expected  for  a  periodic  surface  stress 
(s'-ction  4.1).  In  fact  when  la  =  0(1)  destructive  interference  starts  to  take  place  between  the 
vertical  velocity  induced  by  the  different  patches. 

In  a  subgrid  parameterization  of  the  effect  of  surface  roughness  in  numerical  models,  we  can 
say  that  if  the  rough  patches  are  small  in  comparison  to  the  inverse  of  the  Scorer  parameter, 
the  perturbations  interfere  distructively  and  therefore  need  not  be  explicitly  resolved,  but  can  be 
averaged  beforehand.  If  the  horizontal  scale  of  the  rough  patches  is  equal  or  larger  than  tiie  ir.verse 
the  ambient  Scorer  parameter,  the  vertical  velocity  perturbation  needs  to  be  resolved  explicitly 
or  parameterized  through  the  the  use  of  the  theory  presented  in  this  paper. 
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6  Conclusions 


VVe  have  shown  that  a  horizontal  change  in  surface  roughness  can  induce  substantial  vertical  velocity. 
The  vertical  velocity  can  be  in  the  form  of  a  propagating  waves,  which  can  penetrate  deeply  into 
the  atmosphere.  The  vertical  velocity  can  be  in  the  form  of  trapped  waves  as  well.  In  this  case 
the  perturbation  can  be  physically  relevant,  since  the  maximum  is  placed  at  the  top  of  the  stress 
layer,  i.e.  in  the  region  where  it  is  inportant  to  have  positive  vertical  velocities  in  order  to  trigger 
cumulus  convection.  The  nature  of  this  perturbations  depends  on  the  enviromental  parameter  and 
the  horizontal  destribution  of  the  surface  roughness.  In  general  the  vertical  wave  number  is  closely 
related  to  the  ambient  Scorer  parameter.  The  horizontal  scale  of  the  perturbation  equals  the  scale 
of  the  surface  roughness  when  the  transition  is  smooth.  For  abrupt  changes  in  surface  roughness 
the  horizontal  scale  of  the  perturbation  equals  the  inverse  of  the  ambient  Scorer  parameter. 

When  the  horizontal  distribution  of  the  rough  patches  is  periodic,  resonance  amplification  can 
occur  when  their  wave  number  approaches  the  value  of  the  ambient  Scorer  parameter. 
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7  Appendix 


Tlie  forced  solution  of  eq.  (3.1)  is: 

ko 


^'JoTced{l,z)  -  Jq 


2  I  /  r  ,  ,s\rM Hoik)  (z  -  z'))cos{kx) 


(  G^{k)fMk,z') 

V 


ko  0 


{z-z'))co^kx)GHk)f,Ak,z') 
Po{k)  pU 


OO  2 

jikj 

I  0 


dz 


fS'mhipoik)  (z  -  z’))cos(kx) 

Poik) 


(7.1) 
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Hovever  the  first  and  the  third  term  do  not  satisfy  to  the  trapped  wave  condition,  and  the  second 
term  does  not  satisfy  to  the  radiation  condition  (Sommerfelf,  1912-48).  In  order  to  have  meaningful 
solutions,  which  satisfy  to  the  physical  and  to  the  boundary  constraints,  we  have  to  add  a  number 
of  free  modes. 

Using  Green  function  theory,  we  seek  for  the  solutions,  go{k,z  -  z')  and  giik,z  -  z'},  to  the 
governing  equation  (3.1)  for  a  point  source  forcing  ^(i', z').  If  go{k,z-z')  and  gi(k,z-z')  satisfy  to 
the  boundary  condition,  and  go{k,  z  -  z')  satisfies  to  the  trapped  wave  condition,  while  gi(k,  z  -  z') 
satisfies  to  the  radiation  condition,  the  total  solution  for  a  given  forcing  is; 


foTcedi^ ^ 


G^(k)f,,(k,z') 

pU 


-f 


/"/ 


dz'g,{k,z-z') 


G'^{k)f„{k,z') 

pU 


OO  Z  t 

+  Jdk  Jdz'goik,z-  z')l- 


G^(k)f„{k. 
pU 


I-  rM 


(7.2) 


7.1  Propagating  waves  and  radiation  condition 


The  Green  function  (the  response  to  a  point  source  S{x',z'))  for  the  upward  propagating  wave, 
which  satisfy  to  the  radiation  condition,  is: 

-  ^,^siu{pi(k)(z  -  z'))cos(kx)  +  cosipiik)  {z  -  z'))s\r){kx) 

gup[f^,  <  -  z  )  =  -h  e(z  -  z  ) - - - — -  ,.3) 

X  p^(k) 


_  ,  .  _  ^,.sin(/i|(^)(z- z^)-H;3;) 

X  ‘  /i,(fc) 

The  second  term  in  eq.  (7.3)  is  the  added  free  mode. 

Remark  For  verification,  we  derive  the  boundary  value  Green  Junction,  gg,.{k,z): 

9Bcik-:)  =  -  iim  (':^5up(A’,z  -  r'))  =  -  cos{pi{k}z  kx) 

2  -*0  \  oz  /  j>_o  X 

which  is  the  Green  Junction  Jor  radiative  wave  in  the  mountain  problem. 

1  he  downward  propagating  wave  is: 


9d  ^u.m(k.  z  -  z')  =  -[/fc(z)  -  ne{z  -  z')] 

TT 


sinjpiik)  (z'  -  z)  -I-  kx) 
P\{k) 


(7.4) 


(7.5) 
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and  the  wave  reflected  form  the  ground  is: 

-  ,,  2  „  ,  ,sin(^i(A:)  (z  +  z')  +  A:x) 

g„j{k,z-  z  )  =  -He{z) - — — - 

IT  A'll*) 

The  resulting  Green  function  in  the  propagating  wave  region  is: 

gi{k,z  -  z)  =  gup{k,z  -  z')  +  gdoum{k,z  -  z')  -  grej{k,z  -  z') 


7.2  TTapped  waves 


The  Green  function  for  the  upward  trapped  wave  is: 

.  ,^s\nh(pio{k)  (z  -  z'))  -  cosh{fio(k)  {z - 

gxtp(k,.i.  z  )  —  Tie(z  z  ) 

TT  fi^{k) 


cos{kx] 


tr  t^o{k) 

The  second  term  in  eq.  (7.8)  is  the  added  free  mode. 

Remark  For  verification,  we  derive  the  boundary  value  Green  function,  gg^{k,z): 

d  2 

9bc(^'')  =  -  I™  ^,9up{k,z-  z')  =  -exp(-/io(fc)  (z  -  z'))cos(kx) 

i'— 0  OZ  TT 

u  liich  in  the  Green  function  for  the  trapped  U'ave  in  the  mountain  problem. 

The  downward  trapped  wave  is: 

^  ,^,exp(-/ij^)  (r' -  z))  ,,  , 

9dou-n{k,z  -  z  )  =  --  //c{z)  -  //e(z  -  z  )  - — - cos(kx) 

IT  pjk) 

and  the  wave  reflected  form  the  ground  is: 

-  ,,  M  2  „  ,  ,exp(-/io(^)  (i- +  ^'))  ,,  , 

gref{k,z  -  z  )  =  --He{z) - — - cos(fcx) 

TT  /io(fc) 

The  resulting  Green  function  in  the  trapped  wave  region  is: 


(7.10) 


go{k,Z  -  z')  =  gupik.Z  -  z')  +  gdou:n{k,Z  -  z')  -  grej{k,Z  -  z') 


(7.12) 
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List  of  Figures 


Figure  1.  Contours  of  vertical  velocity  induced  by  a  sudden  change  of  surface  roughness  as  simu¬ 
lated  by  a  non-linear  numerical  model. 

{U  =  0  m  s~^  and  /  =  10“^  m“*) 

Figure  2.  Contours  of  vertical  velocity  induced  by  a  sudden  change  of  surface  roughness  as  simu¬ 
lated  by  a  non-linear  numerical  model, 
f  r  =  10  m  and  /  =  10~^  m~^) 

Figure  3a.  Vertical  velocity  induced  by  a  bell  shaped  surface  stress  in  a  stratified  atmosphere. 

(f.  =  8  Am,  C  =  10  m  and  /  =  10"^ 

Figure  3b.  Approximate  solution. 

{(I  =  8  Am.  U  =  10  m  and  I  =  10“^  m“*) 

Figure  3c.  Approximate  solution  in  a  non-stratified  atmosphere. 

((7  =  8  Am,  U  =  10  m  and  /  =  0.0  m“*) 

Figure  -la.  Vertical  velocity  induced  by  a  bell  shaped  surface  stress  in  a  stratified  atmosphere. 

( (i  =  1  Am,  V  =  10  m  and  /  =  10“^  nr"*) 

figure  ‘lb.  .Approximate  solution. 

la  =  1  Am.  U  =  10  m  and  I  =  10“^  m~*) 

f  igure  5a.  Vertical  velocity  induced  by  sudden  change  of  surface  stress  in  a  stratified  atmosphere. 
[V  =  10  m  5”*  and  /  =  10"^  m~*) 

Figure  5b.  Approximate  solution  in  a  stratified  atmosphere. 

(f  =  10  m  and  I  =  10“^  m~*) 

Figure  5c.  Approximate  solution  in  a  non-stratified  atmosphere. 

[V  =  10  m  and  I  =  0.0  m"*) 

Figure  6.  Waves  in  a  stratified  free  atmosphere,  generated  in  a  non-stratified  stress  layer. 
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{U  =  \0  m  s  =  10"^  m  '  and  I'  =  0.0  m  ') 


Figure  7.  Waves  ducted  within  a  stratified  stress  layer,  capped  by  a  weakly  stratified  free  atmo¬ 
sphere. 

(U  =  10  m  I  =  0.0  and  I'  =  10“^  m~^) 


Figure  8.  Vertical  velocity  induced  by  sudden  change  of  surface  stress  followed  by  smooth  patches 
alternated  with  rough  patches. 

(a  =  0.5, 1,2, 4, 8  km,  U  =  10  m  and  /  =  10~^ 
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pg....  and  ..  refere  to  the  integrals  tables  by  Gradshteyn  and  Ryzhik 


8  Derivation  of  the  Approximations 


Bell  shape  stress 


(a)  -  When  kga  s;  1  or  larger,  /q,  >>  /j  +  /^j,  and  the  wave  is  trapped 

The  contributions  of  the  second  and  third  integrals  are  negligible  in  comparison  to  the  contri¬ 
bution  of  the  first  integral  contribution: 

w{x,z)s:w  /  dkGoik)w^„{k,z)  a  exp(-ak) 

Jo 


-  a  i  r,  i  2ia  +  z-h)^  1 

=  w—  <  He(z  -  h) - :r  -  - - 

ko\  '[^a  +  Z-/!)2  +  z2j2  [(fl  +  2-, 


h)2  +  x'‘ 


+ 


=  I  -  k) 


(fl  +  z  -  ) 

[(a+Z-/l)2  +  l2]2j 


with  Po  =  k  and  Goik)  =  7“  + 


W  hich  is  a  trapped  wave. 


(b)  -  II  hin  la  >>  1  (but  kgO  <  1),  Ii  >>  lo^  +  /oj-  and  the  wave  is  hydrostatic 

The  third  integral  does  not  contributes,  because  the  exponential  decay  of  the  Fourier  transform 
of  the  bell . 

f°° 

U’(i,z)  dkGi{k)w^^{k,z)  a  exp(-aA:)  =  /j 

Jo 

/•'» 

Ii  =  w  a  He{z  -h)  dks]i\{X(z  -  h)  +  kx)exp{-ak)  +  ■  •  • 

Jo 


=  u: 


[a2  +  i2 


\He{z  -  /i)  [asin(A(z  -  h))  +  X  cos(A(z  -  h))]  + 


with  Pi  =  /;  Gi{k)  =  1  and  A  =  \fl^  —  kl 
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(c)  -  ir/jen  la  =  0(1),  <<  Ii  +  loj,  and  the  wave  is  non-hydrostatic 

U'(x,  r)  ss /i  + /oj  =u}|jf  dkGl{k)w^^^{k,z)  a  exp{-ka) -{■  dkGo{k)w^,^{k,  z)  a  exp(-A:a)| 


I\  —  w  j  dkGi{k)v:^,^{k,  z)  a  exp(-ka)  s:  w  a  exp(-la)  f  dkG\{k)w^^^(k,  z) 
Jo  Jo 


u- a  exp{-la)ne(z  -  h)  sin(/( 


-  h))  f 

Jo 


' ^^cos(<:(t  -  z  -ir  h)) 


/p  - 


=  u-  la  exp(-la)He(z  -  h)'^  (sin(/(3  -  h))  J^(/(x  -  x  -  L.)  +  cos{l(z  -  h))  lJg(l[x  -  z+  h))]  + 


with  =  I  -  k  and  G'i(it)  = 


/P  -  ’ 


dkGo(k)w^^{k,z)  a  exp{-ka) 


in  tiic  near  field,  large  wave  numbers....  /q,  can  be  written  in  terms  of  incomplete  F  functions: 
/o.  =  w  la  |/?e(x  -  /i)  |-r(0,/(a  +  ;-/!  +  ix))  +  ^r(0,/(a  +  z  -  h  -  ix))  +  ■  •  • 

with  Pa  =  fc;  and  Go{k)  =  p 

k 

4S9  ^  3.941.4 
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which,  through  the  use  of  the  average  value  theorem,  lo^  can  be  written  as.... 

JfOO 

'  dkexp(—k(a  +  z-h))cos{kx)  +  --- 
0 


=  wHe{z  -  h) 


a{a  +  z  —  h) 
[(a  +  2  -  /i)2  + 


+ 


with  Ho  =  k  and  Go{k)  =  1 


in  the  far  field,  wave  numbers  close  to  I 


Io2  u!  a  exp(-/(a  +  z  -  h)) 


°°  cosffci)  , 
dk—p====  + 


Vk^  -  P 


=  —w  la  He(z  -  h)  exp(-/(fl  +  z  -  h))  +  ■■■ 


Ho  =  /  and  Go(k)  =  , 

^/PTTTI 
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(dj  -  11  heri  la  <<  1  (but  exp{  -koa)  ~  1),  lo^  >>  hi  +  h’  the  flow  is  irrotational 
We  assume  /  =  0.0.  then.... 


I'-i  I 


r^o  roc 

.  2 )  =  7]  -f  7^2  =  u-  /  dkGl(k)w^^^{k.z)  a  exp{—ak)  +  w  /  (7A-G'o(A')u-u„(A:. 2)  a  exp(-flA  ) 

Jo  Jkc 


.Most  of  the  contribution  of  7i  comes  for  k  ss  ko,  but  this  contribution  is  canceled  by  very  rapid 
oscillations,  7]  s:  0.0. 

The  far  field  the  contribution  of  the  trapped  wave  Ioq  comes  for  k  ss  ko,  but  this  contribution 
i.s  canceled  by  very  strong  trapping.  Therefore . 


102  -  u.  I  dkGo{k)iL\,Jk,z)  a  exp{-ak) 


=  ^2  {//.(.- 


ki>ko 


^r(0,(a  +  z  ~  h  +  ir)/ri)  +  ^r(0,(a  +  z  -  h  -  ix)ki) 


+ 
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^,2  1 

with  tio  =  wGo{k)  =  and  =  \f{‘^)ko 
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The  Fs  are  incomplete  F  functions.  However,  through  the  use  of  the  average  value  theorem  and 
since  the  vertical  scale  of  the  perturbation  equals  the  depth  of  stress  layer  and  because  the  near 

field  the  contribution  of  the  trapped  wave  comes  for  k  >>  ko,  then . 

u-' 

/oj  =  -TT-  a  He{z  —  h)  I  dk cos(kx)exp{-{a  +  z  -  h)k)  +  •  •  • 

U  Jo 
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Step  function  stress 


(a)  IFai'c  in  a  stmtified  atmosphere,  /o,  <<  /j  +  lo 


w(x,  z)  SS  Jj  +  Ig 


=  2£  f  dkF[ik)G^(k)  r  dkF'g{k)Go{k) 

jr  (^70 


/,  =  ^-He{z-h)  s\n(Hz-h))  f  dk 
2  TT  Jo 


sin(^(x  -  z  +  h)) 

y/F^ 


-  cos(/(2  -  h))  f  dk 
Jo 


cos(fc(x  -  z  +  h)) 
y/F  -  k^ 


1\  =  —He{z  -  h)  [sin(/(2  -  h))  Hoilix  +  a  -  z  +  h))  -  cos{l{z  -  h))  Jo{l(x  +  a  -  x  +  h))]  + 


Pi  ■=■  I  —  k  and 


Gi(k) 


T2  -  k^ 


in  the  near  field,  large  wave  numbers...  can  be  written  in  terms  of  F  functions.... 

4  =  ^  |//C(2  -  h)  '-1({Z  -  h)  +  1(X  +  h))F(-l,/((2  -h)  +  i(x  +  h))) 
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-  -  /»)  -  i(2  +  /l)]r(-l,/((2  -  /l)  +  i(x  +  A)) 


+ 


U' 


/o,  =  -He{z  -  h) 
TT 


{  [^/((^  -h)  +  t{x  +  h))T{-lM^  -h)  +  iix  +  h))) 


■^![{z  -h)-  i{x  +  h)]r(-l,/((z  -  A)  +  t'(x  +  /2))J  I  + 

with  Mo  =  and  Go(k)  =  ^ 
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u  Jiich  can  be  approximate  to. 


r  rr  ,  ,  ^  _1  3^  +  O 

/oj  =  — //c(z  -  A)tan  ’ - -  + 

X  z  h 


Mo-k  and  Go{k)  =  1 
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in  tiic  far  field,  wave  numbers  close  to  I 
u-2 


Io2  =  ~  J  <^kexp{-k(z  -  A)) 


sin(A(x  +  c)) 
v/F-  /2 


+ 


sf^n(x  +  a)  yo(^(3' +  - /»)exp(/(a  -  A})  • 


and  Mo  =  I  and 


Go{k)  I 


k 
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(b)  Trapped  wave  in  a  non-stratified  atmosphere,  lo^  >>  lo^  +  I\ 


w{x,z)  =  —  r dkF',{k)Go{k) 
IT  Jo 


°°  ,,  s\n(k(x  +  a))  ,  ,  ,  ... 

dk — - ^exp(-ib(2  -  h))  + 

k 


=  ^^He{2  -  /i)tan  *  + 

■kU  z  —  h 


with  po  —  k  and  wGo{k)  =  -jj- 
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Waves  in  a  two  layers  atmosphere 

(n)  Have  in  the  free  atmosphere,  generated  in  a  non-stratified  stress  layer 
.Matching  of  the  solution  emerging  from  the  stress  layer  into  the  free  atmosphere  we  have  a 
trnsniitied  mode  and  a  trapped  mode: 

w{x,  r)  w  /i  +  /oj 

I\  =  —r:He{z  —  h)  I  eos{p\(z  —  h)  kx)  [1  -  exp{-2hk)]  +  •  •  • 

x  t  Jo  k 


t/*  fi\i\(  lcn.\ 

hi  =  —yrHeiz  -  h)  I  dk — - — ^ cos(fci)exp(po(z  -  h))[\  -  exp(-2/iA')]  + 
“kU  Ji  k 

1  c.xpand  the  exponential  in  the  reflected  mode  (emerging  from  the  stress  layer) . 


[1  -  exp(-2/!it)]  =  [l  -  1  +  2/it  -  2{hky  +  ■■■]=  2hk[l  -  hk +  ■  ■  ■] 
furthermore . we  have: 

sin(A:a)  sin( k( a x  -  z  +  h)  +  Hz  -  h))  -  s\n{k{a  -  x  +  z  -  h)  -  l{z  -  h)) 

— - —  cos(/i,(z  -h)  +  kx)  =  — i - — - 
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and 


cos{kx)  =  ■^[sia(k(a  +  i))  -  sin(t(a  -  i)) 

in  the  two  above  eq.s  the  second  term  is  the  downstream  mode,  which  we  neglect.,  then  the 
trasmitted  mode  at  the  upstream  transition  is: 

/i  =  -77  h  ne{z  —  h)  f  dk{\  -  hk) s\n{k(a  +  x  -  z  +  h)  +  l{z  —  h))  +  ■  ■  ■ 

TlJ  Jo 

then  the  trapped  mode  at  the  upstream  transition  is: 

/oj  =  —77  h  He{z  —  h)  f  dk{l  —  hk)sin{k{a  +  x)exp{k{z  —  h))  +  •  •  • 
ttU  Jo 


h  +  lo,  =  -  h) 


^^h[sin(i(j  +  g))  —  s\n{l{z  -  h))]  -  {x  +  a  -  z  +  h)[(l  +  hI)cos(l{x  +  a))  -  cos(Hz  —  h))] 

(x  +  a  -  z  +  /i)2 


+h(x  +  a) 


(x  +  c)^  +  (z  -  h)^  -  2hiz  -  h)' 
[(x  +  a)2  +  (z-h)2f 
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(b)  U’fli'f  ducted  within  the  surface  layer 

1*'  -  The  trapped  mode  in  the  stress  layer  is  trapped  in  the  free  atmosphere  as  well. 

3’’'^  -  The  modes  ducted  within  the  stress  layer  are  shifted  by  2h  at  each  reflection:  the  mode 
refit  cted  by  the  ground  is  shifted  by  2h  at  each  reflection,  the  mode  returning  from  the  top  of  the 
stress  layer  is  shifted  by  2h  at  each  reflection . 

[H,  {=)  -  He{z  -  h)]IZ^,  =  -  He[z  -  h)] 

L  fi 
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Jo  y/l'i  -  fc2 

-  cos(J((2m  +  l)h  -  2))  /'.,.cosW^- ((2m +!)«-»)))' 

Jo  y/fi  -  k} 

[Hc{z)  -  He{z  -  /i)]Q  =  -|-[7/e(z)  -  He{z  -  h)] 

Z  Tt 

s\R(l({2m+l)h+  z))  /'^frSin(<:(T-((2m+  l)/»  +  2))) 

Jo  y/l^  - 

-  cos(/((2m  +  l)ft  +  ,))  + 

Jo  y/l^  -  fc2 

m  =  0,1,2, 3--  =  and  = —2= 

2"=^  ■  The  propagating  mode  in  the  stress  layer  (3”'^),  computed  at  z=h,  is  matched  with  trapped 
in  the  free  atmosphere: 

Matching  of  the  modes  emerging  from  the  stress  layer  with  the  trapped  mode  in  the  free 
it  nuispliere: 

He{z  -  h)[rXg^(z  =z  h}  +  I^jlz  =  h)]  exp(~l{z  -  h)) 

m  =  0,l,2,3  •••  =  / 
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